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Application  of  Satellite  Accelerometer 
Data  to  Improve  Density  Models 


I.  INTRODUCTION 


AFGL  Is  obtaining  comprehensive  solar-cycle  coverage  of  atmospheric  density 

variations  in  the  lower  thermosphere  (between  about  140  and  240  km)  from  satellite 

accelerometer  measurements.  We  accumulated  extensive  low  solar  flux  data  during 

1974  -  1970  from  the  S3-1  and  Atmosphere  Explorer  -C,  -D  and  -E  (AE-C,  -D,  -E) 

experiments.  1  This  data  set  has  been  incorporated  into  an  empirical  model  of 

o 

density  variability  during  low  solar  flux  conditions.  Results  for  moderate  to  high 

3  4 

solur  activity  periods  are  becoming  available  from  the  ROCA  and  SETA  (SETA-1 
and  -2)  Instruments.  Although  only  samples  of  these  data  are  available  for  analysis, 
our  accelerometer  density  data  base  is  being  continually  expanded  as  new  results 


(Received  for  publication  8  August  1984) 
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become  available.  This  data  base  permits  systematic  analyses  of  atmospheric 
density  variability  and  eventual  development  of  improved  models. 

A  statistical  evaluation  of  several  atmospheric  models  using  selected  data  from 
SETA-1,  SETA-2,  and  ROCA  has  been  compared  with  similar  analyses  using  the 
low  solar  flux  data  from  S3-}  and  AE-C,  -D,  and  -E.  Our  present  knowledge  of 
absolute  density  and  its  variations,  based  mainly  on  these  results,  is  discussed 
along  with  a  brief  review  of  the  nature  of  the  problem  and  requirements  for  achieving 
improved  accuracies.  This  preliminary  assessment  is  in  response  to  a  Space 
Division  need  for  a  summary  of  our  present  and  future  capability  to  specify  atmos¬ 
pheric  density  in  the  lower  thermosphere. 

2.  DATA  DESCRIPTION 

Table  1  lists  the  sources  of  the  present  data  set  and  the  dates  of  data  acquisition. 

The  single  axis  MESA  (Miniature  Electrostatic  Accelerometer),  ^  flown  on  AE-C, 

i  o 

-D  and  -E  (triaxial  configuration!  and  on  S3-1,  has  been  described  previously. 

On  the  S3 -4  satellite  the  MESA  unit  was  mounted  on  a  table  which  could  be  rotated 

3 

perpendicular  to  the  satellite  flight  direction.  This  ROCA  (Rotatable  Calibration 
Accelerometer'  experiment  provided  an  orbital  bias  calibration  capability.  The 
SETA  (Satellite  Electrostatic  Triaxial  Accelerometer)  instruments4  permit  simul¬ 
taneous  measurments  of  three  orthogonal  acceleration  components  with  a  single 
proof  mass.  Horizontal  and  vertical  winds,  as  well  as  density,  can  be  determined 
using  only  one  instrument. 


Table  1.  Satellite  Accelerometer  Data  Sources 


Satellite 

Data  Acquisition  Period 

AE-C 

Jan  -  Dec  74 

S3- 1 

Oct  74  -  May  75 

AE-D 

Oct  75  -  Jan  76 

AE-E 

Nov  75  -  Nov  76 

S3-4 

May  -  Aug  7  8 

SETA-1 

Mar  -  Apr  79 

SETA-2 

May  -  Nov  82 

5.  Champion,  K.  S.  W.  ,  and  Marcos,  E.  A.  (1973)  The  triaxial  accelerometer 
system  on  Atmosphere  Explorer,  Radio  Sci.  8:197. 


The  data  sets  of  Table  1  have  been  divided  into  two  groups.  Group  I  consists 

of  the  low  solar  flux  data  from  S3-1  and  AE-C,  -D  and  -E.  Group  II  consists  of 

the  moderate  to  high  solar  flux  data  from  S3-4  (ROCA)  and  SETA-1,  -2.  The 

1  2 

major  features  of  the  Group  I  data  distribution  '  summarized  in  Figures  la-le 
are: 

(a)  Altitude  coverage  135-240  km. 

(b)  Low  solar  flux  conditions,  with  96  percent  of  the  data  at 
Fjq  ^  =  100  units. 

(c)  An  extensive  range  of  geomagnetic  conditions. 

(d)  Good  local  time  coverage,  particularly  from  AE-C  and  AE-E. 

(e)  Latitude  coverage  mainly  at  low  and  middle  latitudes  with  79  percent  of  the 
data  between  the  equator  and  ±  60°  (geographic).  There  is  a  non-uniform  distribu¬ 
tion  of  the  latitudinal  coverage.  High  latitude  data  (>  60°)  occur  mainly  during 
winter  and  spring  over  the  northern  hemisphere  and  during  summer  over  the 
southern  hemisphere. 


ALTITUDE  t KM  3 


S0LRR  FLUX 


Figure  la.  Histogram  of  Group  I  Data 
Distribution  as  a  Function  of  Altitude 


Figure  lb.  Histogram  of  Group  1  Data 
Distribution  as  a  Function  of  Solar  Flux 
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Figure  It;.  illBtogrum  of  Group  I  Data  Figure  Id.  Histogram  of  Group  I  Data 

Distribution  as  u  Function  of  Geomag-  Distribution  as  a  Function  of  Local  Time 

nctlc  Activity 


DAY  OF  YEAR 


Figure  le.  Distribution  of  Group  I  Data  at  200  km  In  Ooographlc 
Latitude  vs  Day  of  Year  Coordinates 
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Tho  available  Group  It  data,  obtained  between  170-240  km,  complement  and 
extend  the  Group  1  data,  Figure  2  shows  that  the  aolar  flux  condition*  extend  over 
u  range  of  almost  200  unite  to  a  maximum  value  of  270  unite.  Therefore,  the  data 
are  for  moderate  to  high  solar  activity.  Figure  3  la  a  hletogram  of  the  diatrlbutlon 
of  geomagnetic  activity  (Kp)  showing  that  a  wide  rango  of  conditions  was  encountered, 
very  much  like  that  for  Group  I  (Figure  lc'.  Figure  4  is  presented  In  order  to 
summarize  the  dramatic  variations  In  solar  and  geomagnetic  condition*  encountered 
during  tho  SETA-2  flight.  The  Fjq,  7  cm  *olar  flux  ar.d  dally  eum  of  the  3-hourly 
Kp  indox  are  shown  as  function  of  day  of  year,  from  l  May  -  7  October  1082,  Dis¬ 
tribution  of  the  data  in  day  of  year  vs  geographic  latitude  Is  shown  In  Figure  5. 
Results  are  obtained  from  the  northern  hemisphere  polar  region  to  about  40*  south 
latitude  from  March  to  November.  The  major  feature  of  the  Group  11  data  arei 

(a)  Altitude  coverage  approximately  170-240  km. 

(b)  Moderate  to  high  solar  flux  conditions, 

(c)  Local  times  mainly  near  1030  and  2230  hours. 

(d)  Nearly  uniform  latitudinal  coverage  from  north  polar  regions  to  about 
40°S  for  the  period  March  -  November, 

The  combination  of  an  extensive  lower  thermosphere  densitv  data  set  obtained 
over  a  wide  latitude  region  with  routine  monitoring  at  high  latitudes  has  not  been 
previously  achieved.  The  present  and  planned  Group  II  data  provide  a  unique  oppor¬ 
tunity  to  Investigate  the  role  of  the  various  heat  sources  that  affect  neutral  atmos¬ 
pheric  density.  For  example,  effects  of  geomagnetic  storms  can  now  be  studied  to 
determine  their  dependence  on  latitude,  season  and  solar  flux  at  morning  and  evening 
local  times.  Further,  any  proposed  Improved  indicators  of  solar  and  geomagnetic 
activity  can  be  readily  evaluated, 


Mr— — r— — •'  r . 1  -  **r*  1  —r  1  1 1  —  i  - 1 


SOLAS  FLUX  (FI0.7J 


Figure  2,  Histogram  of  Group  II  Data  Distribution  as  a  Function 
of  Solar  Flux 
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Figure  4.  Solar  Flux  and  Geomagnetic  Activity  Indicators 
vs  Day  of  Year  for  SETA-2  Data 
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MY  OF  YEAR 


Figure  5.  Distribution  of  Group  II  Data  in  Geographic  Latitude  vs  Day  of 
Year  Coordinates  .  , 


Figure  6  outlines  the  data  selection  and  analysis  performed  for  this  study.  The 
satellites,  time  periods  of  measurements  and  total  number  of  orbits  for  which  data 
are  available  are  given  in  the  top  boxes.  Portions  of  each  data  set  utilized  are 
shown  in  the  next  row.  Various  sub-sets  were  studied  and  are  specified  as  required 
in  the  next  section.  Software  was  developed  to  calculate  ratios  of  measured  density 
to  numerous  models  in  the  "Model  Ratio  Program."  The  models  considered  were 
six  versions  by  Jacchia:  J64,  J70,  J71,  J73,  and  JWB  (see  Reference  6  and  refer- 
ences  therein)  and  J77,  the  MSIS  (Mass  Spectrometer  and  Incoherent  Scatter)  model 

without  (MSIS77'**  and  with  (MSIS79'^  longitude  variations,  the  U.  S.  Standard  Atmos- 

6  6 
phere  1962,  the  U.  S.  Standard  Atmosphere  Supplements,  1966,  the  Lockheed-NASA 

C  C 

model  and  DENSEL.  The  ratios  were  used  as  input  to  a  "Statistics  Program"  for 
-  calculations  of  mean  values  and  standard  deviations  in  selected  Kp,  latitude  and 
local  time  bins. 
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J,  Geophys.  Rea.  JJj2:2139. 

9.  Hedin,  A.  E. ,  Reber,  C.A.,  Spencer,  N.W.,  and  Brinton,  H.  C.  (1979' Global 

model  of  longitude/UT  variations  in  thermospheric  composition  and 
temperature  based  on  mass  spectrometer  data,  J,  Geophys.  Res.  84:1. 
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Figure  6.  Satellite  Accelerometer  Density  Data  Sources  Utilized 
for  Statistical  Analyses 


Digital  filters  were  used  to  separate  drag  from  extraneous  noise  accelerations 
in  deriving  atmospheric  density  from  accelerometer  data.  These  filters  differed, 
depending  mainly  upon  the  spacecraft  spin  rate.  Comparison  of  standard  deviations 
between  different  data  sets  requires  consideration  of  the  data  smoothing  due  to  the 
filtering.  For  the  Atmosphere  Explorer-C,  -D  and  -E  satellites,  data  obtained  in 
the  "despun"  mode  (one  revolution  per  orbit)  were  reduced  using  a  "l0-15"-sec  filter. 
This  filter  rejects  all  data  with  a  period  less  than  or  equal  to  10  sec  and  passes  all 
data  with  a  period  equal  to  or  greater  than  15  seconds.  Between  these  two  values 
the  filter  "transmission  coefficient"  increases  from  zero  to  100  percent  with  a 
quarter  sine  wave  shape.  Data  obtained  in  the  "spinning"  mode  (nominally  one  revo¬ 
lution  per  15  sec)  were  reduced  using  a  "30-77  "-sec  filter.  Therefore,  assuming 
a  satellite  velocity  of  8  km/ sec,  atmospheric  phenomena  with  wavelengths  less  than 
about  80  km  (despun  data)  or  240  km  (spinning  data)  are  smoothed  out  by  the  digital 
filtering.  Data  from  S3-1  were  processed  by  first  applying  a  "notch"  filter  to  re¬ 
move  noise  with  periods  between  13  and  50  sec,  calculating  density  as  a  function 
of  time  and  then  applying  a  "4-8"-sec  filter.  For  the  Group  II  data  a 


10.  Noonan,  J.  P. ,  Fiorettl,  R.  W. ,  and  Hass,  B.  (1975)  Digital  Filtering  Analysis 
Applied  to  the  Atmosphere  Explorer-C  Satellite  MESA  Accelerometer  Data, 
AFC  RL-TR-75- 0293,  ATT  A 0 137 5 5/ .  . 
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";tO- 100" -sec  filter  was  generally  uaed.  A  special  study  of  higher  frequency 
atmospheric  structure  was  also  conducted  UBing  ten  days  of  SETA-2  data  processed 
with  a  "  15-20"-soc  filter. 


3.  RESULTS 

3.1  Group  I  Data  Mean  Ratio*  and  Standard  Deviation* 

Table  2  provides  a  statistical  evaluation  of  the  entire  AE-C,  -D  and  -E  and 
S3-1  data  base.  As  noted  in  Section  2,  the  quantity  statistically  analyzed  Is  the 
ratio  (R)  between  the  measured  density  and  model  density.  Results  are  given  for 
mean  ratios  (IT): 

N 

E  R  (1) 


and  standard  deviations  (ct ^): 


1 

N  -  1 


N 

E 


n=  1 


(R-TV) 


(2) 


where  N  is  the  number  of  data  points.  Models  used  for  this  study  are  listed  in  order 
of  increasing  percent  standard  deviation.  For  the  first  three  models  listed,  the 
mean  ratios  are  within  1.  045  ±  0.  025  and  the  standard  deviations  are  within 
14.  5  ±  0.  3  percent,  indicating  excellent  overall  agreement.  Note  that  the  unpub¬ 
lished  J7  3  model  has  the  lowest  standard  deviation,  0.  6  percent  less  than  that  of 
J71.  At  the  other  extreme,  the  root-mean- square  (rms)  variations  for  the  DENSEL 

ratio  are  exceeded  only  by  US62  for  this  low  solar  flux  data  set. 

2 

The  empirical  Atmospheric  Density  Model  (ADM)  developed  from  the  above 
data  set  has  been  used  to  evaluate  other  model  depictions  of  density  variability  with 
solar  flux,  geomagnetic  activity,  local  time  and  day  of  year.  As  an  example  of  that 
analysis.  Figures  7a-7c,  taken  from  Reference  2,  compare  the  local  time  de¬ 
pendence  of  density  from  ADM  with  that  of  the  J71,  J77  and  MSIS77  models.  At 
160  km  the  local  time  behavior  is  dominated  by  the  semidiurnal  tide,  in  very  good 
agreement  with  the  MSIS  model  and  contrary  to  the  J71  and  J77  diurnal  tide  pre¬ 
dictions.  At  200  km  and  240  km  the  diurnal  tide  becomes  predominant,  but  appears 
to  be  overestimated  in  J7t  and  J77  and  underestimated  in  MSIS.  A  statistical 
evaluation  of  ADM,  using  the  total  data  set  it  is  based  upon,  provided  the  following 
result: 
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Model 


Mean  Ratio 


Std.  Dev.  (%> 


ADM  1.000  12.8 

This  standard  deviation  is  1.4  percent  smaller  than  that  of  J73.  The  calculated 
standard  deviation  represents  an  Internal  test  of  the  consistency  between  the  mathe¬ 
matical  representation  and  the  data  set  used  for  model  development.  The  actual 
test  of  a  model  is  Its  ability  to  represent  Independent  data  (not  used  in  its  construc¬ 
tion).  However,  preliminary  tests,  using  a  model  constructed  from  a  portion  of 
the  data  and  evaluating  It  with  the  remainder  of  the  data,  have  indicated  the  Improved 
accuracy  of  ADM.  These  tests  included  evaluating  a  model  based  on  cnly  northern 
hemisphere  data  with  data  obtained  in  the  southern  hemisphere  and  evaluating  a 
model  based  on  three  satellite  data  sets  with  results  from  the  fourth  data  set. 

Future  efforts  directed  toward  incorporating  the  Group  II  data  into  an  Improved 
version  of  the  ADM  model  can  lead  to  increased  understanding  of  density  variability 
at  high  latitudes  and  high  solar  flux  conditions. 


Table  2.  Group  I  (AE-C,  -D,  -E  and  S3-1)  Total 
Mass  Density  Ratio  to  Models  (Altitude  =  135-240  km; 
Local  Time  =  0-24  hr;  Kp  Range  =  0-9) 


Model 

Mean  Ratio 

Std  Dev 
(%) 

1. 

J73 

1.057 

14.  2 

2. 

MS1S  (77' 

1.  021 

14.6 

3. 

J71 

1.  070 

14.  8 

4. 

J77 

1.  032 

15.6 

5. 

J70 

1.  022 

15.  8 

6. 

US66 

0.  931 

16.  1 

7. 

L-NASA 

0.  910 

16.6 

8. 

J64 

0.  921 

17.  3 

9. 

JWB 

1.  000 

19.6 

10. 

DENSEL 

1.  214 

28.  1 

11. 

US62 

0.  773 

33.6 

Total  No.  of  Points  =  154,749 
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Figure  7a.  Density  Variation  at  160  km 
as  a  Function  of  Local  Time  (Day  =81, 
Latitude  =0°,  F  =  Tr  =  75,  Kp  =  0) 


LOCAL  TIME 


Figure  7b.  Density  Variation  at 
200  km  as  a  Function  of  Local  Time 
(Day  =  81,  Latitude  =0°,  F  =  ¥  =  75, 
Kp  =  0) 
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LOCAL  TIME 


Figure  7c.  Density  Variation 
at  240  km  as  a  Function  of  Local 
Time  (Day  =81,  Latitude  =  0°, 

F  =  T  =  75,  Kp  =  0) 


3.  1.  1  FREQUENCY  AND  PROBABILITY  DISTRIBUTIONS 

Histograms  of  the  frequency  distribution  and  curves  of  the  probability  distribu¬ 
tion  of  the  data,  shown  as  percent  departure  from  models  were  given  for  Group  I 
data  in  Reference  1.  That  study  included  calculations  of  the  third  and  fourth 
moments  (skewness  and  kurtosls,  respectively)  as  well  as  the  mean  and  standard 
deviation.  Although  results  showed  some  deviations  from  a  normal  distribution, 
there  are  no  infrequently  occurring  very  large  differences.  Data  for  AE-C  and 
AE-D  are  reproduced  in  Figures  8a  and  8b,  and  their  cumulative  probability 
distributions  are  shown  in  Figures  9a  and  9b.  Similar  analyses  are  planned  using 
the  Group  II  data. 

3.  1.  2  MODIFICATION  OF  GROUP  I  DATA  BASE 

The  Group  I  data  base  was  modified  by  considering  only  Atmosphere  Explorer 
data  and  limiting  results  to  the  region  170-220  km  for  local  times  08-24  hr.  This 
was  done  to  provide  a  data  set  that  was  more  self-consistent  and  more  compatible  for 
later  comparison  with  the  Group  II  data.  The  AE  satellites  had  similar  orbital 
parameters,  shapes,  drag  coefficients  and  error  sources.  Density  values  were 
derived  with  identical  filters.  The  170-220  km  region  is  also  used  for  the  Group  II 
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data  for  this  report  and  all  Group  II  data  points  have  local  times  between  08  and 
24  hr,  although  they  are  heavily  concentrated  near  10  and  22  hrs  (limiting  the 
Atmosphere  Explorer  data  base  to  a  similar  local  time  distribution  would  result 
in  elimination  of  about  90  percent  of  the  available  data).  Statistical  results  given 
in  Table  3  for  this  limited  Group  I  data  set  show  MSIS77  has  the  smallest  standard 
deviation,  1  percent  less  than  its  value  for  all  data  (Table  2).  For  the  first  three 
models  of  Table  3  the  standard  deviations  are  reduced  to  13.9  ±0.3  percent.  Thus, 
the  restricted  altitude  and  local  time  range  results  in  a  reduction  of  at  least 
0.  8  percent  in  the  standard  deviation. 


Figure  8a.  Frequency  Distribution  of  Percent  Difference  of 
AE-C  Data  From  J71  Model 
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Figure  8b.  Frequency  Distribution  of  Percent  Difference  of 
AE-D  Data  From  J71  Model 
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Figure  9a.  Probability  Distributions  for  AE-C 
Data  Shown  as  Percent  Departure  From  the 
J71  Model 
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Figure  9b.  Probability  Distributions  for  AE-D 
Data  Shown  as  Percent  Departure  From  the 
J71  Model 


Table  3.  Modified  Group  I  (AE-C,  -D  and  -E)  Total 
Mass  Density  Ratio  to  Models  (Altitude  ■  170-220  km; 
Local  Time  =  08-24  hr;  Kp  Range*  0-9) 


Model 

Mean  Ratio 

Std  Dev. 

(%' 

1. 

MS1S  77 

1.  016 

13.6 

2. 

J73 

1.033 

13.8 

3. 

J71 

1.  052 

14.2 

4. 

US66 

0.  892 

15.  2 

5. 

J77 

1.  008 

15.  2 

6. 

J64 

0.  890 

15.3 

7. 

J70 

0,  971 

15.8 

8. 

JWB 

0.  927 

16.  2 

9. 

L-NASA 

0.  899 

16.6 

10. 

use  2 

0.  717 

24.  7 

11. 

DENSEL 

1.  233 

28.  3 

Total  No.  of  Points  =  34,819 
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3.  1.  3  EFFECT  OF  DIGITAL  FILTERING  ON  GROUP  I  DATA 


Wave  structures  are  presently  not  Included  In  thermospheric  models.  Digital 
filtering  tends  to  smooth  the  accelerometer  data  with  a  consequent  underestimation 
of  the  variability.  This  factor  has  been  examined  by  partitioning  the  Atmosphere 
Explorer  results  used  for  Table  4  Into  "spun"  and  "despun"  orbits.  These  were 
obtained  on  an  alternating  schedule  during  the  AE  missions  and  thus  provide  com* 
parable  data  sets  as  a  function  of  flux,  Kp,  day,  latitude,  local  time  and  altitude. 
Results  are  given  In  Table  4  using  the  J71,  J77  and  MSIS77  models.  The  despun 
data  filtering  results  In  about  a  1.  7  percent  Increase  In  the  standard  deviation  com* 
pared  to  the  spun  data  and  a  0,8  percent  Increase  compared  to  all  data.  It  Is 
believed  that  any  wave  structures  not  observed  with  the  "10-15"-sec  filter  provide 
negligible  contribution  to  the  variability.  Therefore,  Tables  2  and  4  Indicate  that 
the  best  models  reflect  density  variability  for  this  data  set  with  a  standard  deviation 
of  about  15  percent. 


Table  4.  Comparison  of  Modified  Group  I  (AE-C,  -D  and  -E)  Spun  and  Despun 
Mode  Total  Mass  Density  Ratio  to  Models  (Altitude  ■  170-220  km; 

Local  Time  =  08-24  hr;  Kp  range  ■  0-9) 


All 

Spun 

DeBpun 

Model 

Mean  Ratio 

Std  Dev 

Mean  Ratio 

Std  Dev 

Mean  Ratio 

Std  Dev 

(%) 

(%) 

(%) 

J71 

mm 

14.  2 

mm 

wmm 

HU 

J77 

— 

15.  2 

m 

■8H 

HI 

MSIS 

1.  018 

13.  6 

1 

12.  6 

1.  030 

14.4 

No.  of  Points:  34,819 


17,725 


17,094 


3.2  Group  II  Data  Mean  Ratioi  and  Standard  Devlationi 

Results  of  a  statistical  evaluation  of  the  SETA-1,  SETA -2,  and  S3 -4  accelerom¬ 
eter  measurements  reduced  with  a  "50-100"-sec  filter  are  shown  in  Table  5.  Listing 
of  the  models  is  in  order  of  increasing  average  (unweighted)  standard  deviation 
for  the  three  sets.  The  smallest  standard  deviation  is  associated  with  the  unpub¬ 
lished  J7  3  model,  and  varies  from  7. 70  to  10.  34  percent.  The  latter  value  may  be 
larger  because  of  the  preliminary  state  of  processing  needed  to  incorporate  83-4 
data  into  this  report.  However,  the  standard  deviations  of  8-9  percent  for  SETA-1 
and  SETA- 2  are  typically  3-4  percent  less  than  the  results  for  the  AE  "spun"  data 
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in  Table  4.  There  are  at  least  two  reasons  why  smaller  standard  deviations  would 
be  expected  from  the  Group  II  data.  First,  the  "50-100"-sec  filter  smoothes  out 
more  of  the  atmospheric  structure  than  the  "30-77" -sec  filter.  For 
the  AE  data  the  difference  between  "spun"  ("30-77"  filter)  and  "despun"  ("10-15" 
filter)  data  was  1.  7  percent.  A  similar  difference  might  be  expected  between  the 
Group  I  and  Group  II  data  because  of  the  additional  range  of  wavelengths  eliminated 
by  the  "50-100"-sec  filter  and  the  increase  in  occurrence  and  amplitude  of  longer 
period  waves.  11  Second,  the  Group  I  data  are  for  a  mix  of  altitude,  latitude  and 
local  time  conditions,  whereas  for  the  Group  II  data,  these  conditions  are  highly 
correlated.  For  example,  an  analysis  of  AE-E  data1  showed  that  the  standard 
deviation  decreased  by  about  2  percent  when  the  data  were  restricted  to  local  times 
between  0800  and  0900  hours. 


Table  5.  Group  II  (SETA-1,  SETA-2  and  S3-4)  Total  Mass  Density  Ratio  to 
Models  (Altitude  =  170-220  km;  Local  Time  =  08-24  hr;  Kp  Range  =  0-9) 


SETA 

-1 

SETA 

-2 

S3  -4 

Model 

Mean  Ratio 

Std  Dev 
(%> 

Mean  Ratio 

Std  Dev 
<%> 

Mean  Ratio 

Std  Dev 
(%) 

1. 

J7  3 

0.  916 

9.  03 

0.  942 

7.  70 

0.  934 

10.  34 

2. 

J7  0 

0.  994 

8.  79 

0.954 

8.  20 

0.  925 

10.  43 

3. 

J7  1 

0.  935 

9.  10 

0.  983 

8.  07 

0.  974 

10.  68 

4. 

MSIS77 

0.  974 

10.  20 

0.  969 

8.  02 

0.  970 

11.30 

5. 

US  6  6 

1.  009 

10.  30 

0.  894 

8.  84 

0.  877 

10.  78 

6. 

J  64 

1. 010 

10.  41 

0.  889 

8.  85 

0.  870 

10.  7  1 

7. 

MSIS79 

0.  973 

10.  51 

0.  969 

8.  32 

0.  968 

1 1.  65 

8. 

JWB 

0.  913 

9.  84 

0.  829 

9.  07 

0.  827 

12.  17 

9. 

J77 

0.  879 

12.  26 

0.  925 

10.  06 

0.  979 

11.  65 

10. 

L-NASA 

0.  986 

9.  81 

0.  914 

9.  36 

0.  893 

12.  17 

11. 

US  6  2 

1.  158 

11.  17 

0.  973 

12.  12 

0.  927 

16.  97 

12. 

DENSER 

0.  890 

14.  32 

0.  938 

11. 44 

0.  770 

17.  74 

No.  of  Points  42.835  32.950  36,953 


11.  Hines,  C.O.  (19741  The  Upper  Atmosphere  in  Motion,  American  Geophysical 
Union  Monograph  18. 
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There  are  gome  differences  in  the  relative  accuracies  of  the  models  with 
respect  to  the  three  data  sets  of  Table  5,  Tho  J7  1  model  provides  a  better  repre¬ 
sentation  of  density  variations  encountered  by  SETA-1  and  S3-4  than  does  MSI877. 

For  the  SETA -2  data,  MSIS77  is  slightly  better  than  J7  1,  The  J77  model  apparently 
is  not  an  improvement  over  the  other  Jacchia  models,  It  shows  a  standard  deviation 
2-3  percent  greater  than  J7  1  and  has  an  overall  ranking  of  fourth  from  last  for 
the  three  data  sets,  It  is  not  recommended  for  operational  applications. 

Mean  ratios  in  Table  5  arc  less  than  one  for  the  most  accurate  models  whereas 
the  same  models  had  mean  ratios  greater  than  one  for  the  Group  I  data.  Since  .the 
drag  coefficients  for  both  the  Group  I  and  Group  II  satellites  are  estimated  to  be 
within  t  10  percent,  14  the  difference  cannot  be  unambiguously  related  to  solar  flux. 
The  mean  values  for  the  first  four  models  listed  in  Table  5  are  0.  058  *  0,  020, 

0,  002  ±  0.  021  and  0.  051  ±  0.  025  for  SETA-1,  SETA-2  and  S3-4,  respectively.  The 
variations  are  similar  to  those  found  for  the  Group  I  data.  Interestingly,  MSIS77 
shows  least  change  in  mean  value  for  the  three  Group  II  data  sets.  More  extensive 
analyses,  as  accomplished  in  Reference  1,  are  planned  for  the  larger  data  sets  to 
be  acquired.  These  studies  will  establish  quantitatively  the  relationship  between 
density  and  solar  flux. 

3.3  Data  Partitioning  Into  Kp  Bini 

A  further  analysis  of  the  J71,  J77  and  MSIS77  models  has  been  carried  out  by 
partitioning  the  data  into  three  Kp  bins;  Kp  values  0  to  3,  3+  to  4+  and  5-  to  9. 

This  study  involved  the  modified  Group  I  data  (Table  4)  and  the  Group  II  data  (Table  5). 

3.  3.  1  GROUP  I  DATA 

Table  0  shows  results  for  the  Group  I  data  (Table  4)  with  spun  and  despun 
orbits  considered  separately.  The  differences  between  the  standard  deviations  of 
the  despun  data  and  the  spun  data  are  larger  for  the  5-  to  9  bin  than  for  the  0  to  3  bin. 
Those  differences  ore  -bout  0,  5  percent,  1.  1  percent  and  0,5  percent  for  the  J71, 
MSIS77  and  J77  models,  respectively.  This  difference  is  presumed  due  to  addi¬ 
tional  wave  structures  associated  with  geomagnetic  activity.  The  J7  1  and  MSIS77 
models  show  little  variation  of  standard  deviation  with  Kp.  MSIS77  actually  shows 
a  slight  decrease  with  increasing  Kp,  but  the  corresponding  values  of  mean  ratios 
show  about  a  10  percent  Increase.  The  indicated  variability  (standard  deviation) 
is  not  representative  of  global  conditions  because  the  data  base  used  for  this  study 
is  weighted  toward  lower  latitudes  where  geomagnetic  storm  effects  are  less 
significant.  The  only  data  poleward  of  07°  geographic  latitude  are  provided  by  tha 
short-lived  AE-D  mission.  Further,  the  limited  data  in  the  high  Kp  bin  are 
representative  of  only  a  few  storms.  Hence  the  altitudes,  latitudes  and  local 
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times  are  more  correlated  than  for  data  tn  the  lower  Kp  bine,  Thl*  lend*  to  pro¬ 
vide  a  standard  deviation  that  la  underestimated  at  high  Kp  bine  relative  to  low  Kp 
bins,  A n a  1  y h  1  n  of  Group  II  data  (Section  3,  3.  2)  eupporte  this  hypotheele, 


Table  6.  Modified  Group  I  (AE-C,  -D  and  -E)  Spun  and  Deepun  Modej  Total 
Muss  Density  Ratio  to  Models  Partitioned  Into  Kp  Bine  (Altitude  ■  170-220  kmj 
Local  Time  ■  08-24  hrj  3  Kp  Bins! 


A,  Model  ■  J7  1 

Spun 

Deepun 

Kp 

Rango 

No. 

Points 

Mean  Ratio 

Std  Dev 
(%) 

No. 

Points 

Mean  Ratio 

Std  Dev 
<%) 

All 

17,725 

1,040 

mam 

KBjZZI 

mm 

0  to  3 

11. 688 

1.044 

EK9 

BB9 

1. 063 

DB 

to  4+ 

4.360 

1.032 

13.66 

1,058 

15.51 

5-  to  0 

1.  877 

1.034 

12.97 

1 H 

1,078 

14.96 

B  Model  •  MSIS77 

Spun 

Despun 

Kp 

Range 

No. 

Points 

Mean  Ratio 

Std  Dev 
(%) 

No. 

Points 

Mean  Ratio 

Std  Dev 
(%> 

All 

17,725 

1.  004 

12.57 

17,094 

1,030 

14,44 

0  to  3 

12,35 

11,  118 

1.002 

14,09 

3»  to  4+ 

8KB 

1.033 

12.42 

4,222 

1,070 

13,67 

■  S -  to  0 

1,  677 

1.077 

10.78 

1,754 

1.  108 

13.64 

C.  Model  -  J77 

Spun 

Deepun 

Kp 

Range 

No. 

Points 

Mean  Ratio 

Std  Dev 
<%> 

No. 

Points 

Mean  Ratio 

Std  Dev 
(%> 

All 

14,32 

17,004 

1.019 

16,01 

0  to  3 

■fHHi 

0.907 

13,66 

11.  118 

1.015 

15.  12  * 

3+  to  4+ 

■hi 

0.  907 

15,  12 

4,222 

1,018 

17.  15 

0-  to  0 

0.  907 

16,32 

1,754 

1,042 

16.  15 

25 


3.3.2  GROUP  II  DATA 


Table  7  shows  results  obtained  by  putting  the  Group  II  data  into  the  same  Kp 
bins.  For  the  three  data  sets  the  increases  In  standard  deviation  from  the  0  to  3 
to  5-  to  9  Kp  bin  range  from  2.  7 -3.  4  percent  for  J7  1,  4.  1-5.5  percent  for  J77  and 
0.  7  to  5.  6  percent  for  MSIS.  Mean  ratios  generally  decrease  for  J7  1  and  J77 
(t'/plcally  3-7  percent),  This  increasing  departure  from  unity,  * 

indicating  that,  on  average,  these  models  overestimate  the  observed  geomagnetic 
response.  MSIS77  shows  a  general  increase  in  mean  ratio  except  for  the  3-*-  to  4+ 
to  5-  to  9  bins  for  the  S3 -4  data.  Results  obtained  from  the  more  uniform  latitudinal 
coverage  of  the  Group  II  data  show  increased  errors  in  the  models  during  geomag- 
netically  disturbed  periods.  However,  the  increases  in  these  errors  may  be  under¬ 
estimated.  Wave  structure  occurrences  are  known  to  increase  during  storm 
periods,  and,  as  noted  previously,  the  relatively  heavy  filtering  of  the  Group  II 
data  means  that  part  of  the  spectrum  of  the  waves  were  not  available  for  the  analysis. 
(The  filter  was  selected  to  permit  convergence  of  an  iteration  technique  to  derive 
density  and  winds  simultaneously  from  SETA  data.  )  Hence,  the  Group  I  data  may 
underestimate  the  standard  deviations  because  of  limited  latitudinal  distributions, 
whereas  the  Group  II  data  underestimate  the  standard  deviations  because  of  data 
smoothing.  A  more  comprehensive  quantitive  understanding  of  atmospheric 
variability  during  low  solar  flux  conditions  will  require  a  more  extensive  set  of 
global  measurements.  The  effects  of  smoothing  on  the  Group  II  data  are  further 
discussed  in  Section  3.  5. 

3.4  Data  Partitioning  in  Kp,  Latitude  and  Local  Time  Bin* 

More  detailed  analyses  of  the  statistical  properties  of  the  accelerometer  data 
have  been  made  by  further  partitioning  data  into  10°  geographic  latitude  bins  and 
day/night  local  time  bins  as  well  as  the  three  Kp  bins  used  previously.  Local  time 
bins  of  08-16  hr  and  16-24  hr  were  used  to  separate  Group  II  dayside  and  nightside 
northern  hemisphere  data  occurring  at  the  same  latitude. 

Mean  values  (relative  to  the  J7  1  model)  as  a  function  of  geographic  latitude  for 
the  SETA-1,  SETA -2  and  S3 -4  data  are  shown  in  Figure  10.  Results  are  given  here 
only  for  data  obtained  during  low  geomagnetic  conditions  (Kp  £  3).  The  satellite 
altitudes  as  a  function  of  latitude  are  similar  for  the  data  sets  plotted.  From  rigfht 
to  left  on  the  graph  the  altitude  decreases  from  approximately  220  km  at  low  northern 
hemisphere  latitudes  on  the  nightside  to  about  170  km  near  30®  N  on  the  dayside.  It 
then  increases  to  about  220  km  at  midlatitudes  in  the  southern  hemisphere.  All 
three  data  sets  show  a  density  bulge  located  around  60-70®  latitude  on  the  dayside. 
This  implies  that  a  localized  high  latitude  density  enhancement  is  a  permanent 
feature  of  the  dayside  auroral  region.  The  causative  mechanism  may  be  heating 
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related  to  the  polar  cusp.  More  detailed  studied  of  this  presently  unmooelled 
tiourcc  of  density  variability  are  being  made.  This  result  does  emphasise  that  the 
combination  of  accelerometer  capabilities  and  orbital  characteristics  provides  a 
unique  opportunity  to  obtain  new  fundamental  information. 


Table  7.  Group  II  (SETA-1,  SETA-2  and  S3-4)  Total  Mass  Density  Ratio  to 
Models  Partitioned  Into  Kp  Bins  (Altitude  ■  170-220  km{  Local  Time  ■  08-24  hr; 
3  Kp  Bins) 


Satellite  ■ 

SETA-1 

Satellite  ■ 

SETA -2 

Satellite  * 

S3 -4 

A. 

Model  « 

J7  1 

Kp 

Range 

Mean  Ratio 

Std  Dev 
<%) 

Mean  Ratio 

Std  Dev 
(%) 

Mean  Ratio 

Std  Dev 
(%) 

All 

9.  10 

0.  960 

7.76 

0.981 

11.05 

0  to  3 

7.  30 

0.971 

7.  14 

1.001 

9.  28 

3+  to  4+ 

0.  925 

7.90 

0.933 

7.49 

1.000 

10.  62 

5-  to  9 

0.  889 

10.  06 

0.  940 

9.45 

0.  932 

12.74 

B. 

Model  » 

MSIS77 

Kp 

Range 

Mean  Ratio 

Std  Dev 
(%> 

Mean  Ratio 

Std  Dev 
<%) 

Mean  Ratio 

Std  Dev 
(%> 

--- 

All 

BB 

10.  20 

ns 

8.02 

11.30 

0  to  3 

9.76 

— 

7.59 

9.22 

3+  to  4+ 

10.  53 

0.990 

7.90 

12.25 

5  -  to  9 

■9 

10.49 

0.990 

9.06 

0.940 

14.82 

C. 

Model  ■ 

J77 

kP 

Range 

Mean  Ratio 

Std  Dev 
(%> 

Mean  Ratio 

Std  Dev 
(%) 

Mean  Ratio 

Std  Dev 

fa) 

All 

ns 

12.  26 

BB 

10.06 

BB 

11.65 

0  to  3 

11.  55 

BB 

8.  62 

8.  97 

3+  to  4+ 

■ns 

12.  89 

0.  908 

11.45 

BB 

11.82 

5-  to  9 

0.  847 

15.  63 

0.875 

12.  62 

0.944 

14.46 
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Figure  10.  Mean  Values  of  Ratios  of  SETA-1,  SETA-2  and  S3 -4  Density 
Data  to  J7  1  Model  Plotted  as  a  Function  of  Geographic  Latitude 
(Kp  S  3) 


Analysis  of  data  obtained  during  low  solar  flux  (Figure  7b)  conditions  showed 
that  the  J7  1  densities  were  underestimated  near  200  km  for  local  times  around 
2200  hr.  This  could  be  the  cause  of  the  observed  mean  density  (ratio)  increase 
with  decreasing  latitude  on  the  nightside  since  the  local  time  variations  are  less 
important  at  high  latitudes.  Further  analyses  of  the  SETA-2  and  S3-4  data  are 
required  to  determine  whether  the  apparent  ratio  increase  with  latitude  in  the 
southern  hemisphere  is  related  to  latitude,  altitude  or  a  systematic  bias  error 
related  to  the  preliminary  reduction  state  of  these  data  sets. 

Latitudinal  changes  of  the  standard  deviations  for  the  Group  II  data  are  shown 
for  all  three  Kp  bins  in  Figures  11a  to  lie.  These  figures  show  that  a  maximum 
error  in  the  models  occurs  near  the  pole,  even  during  geomagnetically  quiet  condi¬ 
tions.  The  high-latitude  standard  deviations  tend  to  increase  with  geomagnetic 
activity.  Differences  in  the  dependence  on  Kp  for  the  three  data  sets  are  con¬ 
sidered  mainly  a  manifestation  of  the  inadequacy  of  Kp  as  an  indicator  of  geomagnetic 
activity.  At  lower  latitudes  the  SETA -2  data  indicate  little  variation  of  the  standard 
deviation  with  Kp,  whereas  the  SETA-1  and  S3-4  results  show  an  increase  in 
standard  deviation  for  the  high  Kp  bin.  However,  analysis  of  the  limited  SETA-2 
data  reveals  that  SETA -2  did  not  encounter  as  large  a  range  of  geomagnetic  condi¬ 
tions  as  SETA-1  and  S3 -4.  Consequently,  a  more  comprehensive  study  of  the 
standard  deviation  vb  latitude  will  be  made  using  smaller  Kp  bins. 
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Figure  11a.  Standard  Deviations  of  Ratios  of  SETA-1  Density  Data  to  J7  1  Model 
Plotted  as  a  Function  of  Geographic  Latitude 

A  similar  analysis  of  standard  deviation  vs  latitude  for  low  solar  flux  condi¬ 
tions  was  made  using  S3-1  data.  This  satellite  was  selected  because  it  had  a  sun- 
synchronous  orbit  and  therefore  provided  good  latitudinal  coverage  with  little  local 
time  variation.  Because  of  polar-crossings,  two  perigee  local  times  ~  22  hr  and 
~08  hr  are  mainly  encountered  at  all  latitudes.  Figure  12  shows  structured  fluctua¬ 
tions  in  standard  deviation  that  are  apparently  related  to  unmodelled  density 
structures.  In  general,  at  high  latitudes  a  is  greatest  for  the  high  Kp  case.  Lowest 
values,  8.  5  ±  0.  5  percent,  occur  between  10°  N  and  30°  N.  From  30-60°  N,  the" 
standard  deviation  is  roughly  11  ±  1.  5  percent  and  from  30-60°  S  it  is  about 
13  ±  3  percent.  The  lack  of  a  more  cohesive  latitudinal  structure  for  indicated 
variability  is  attributed  to  the  mixing  of  altitude,  latitude,  day  of  year  and  local 
time. 
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Figure  lib.  Standard  Deviations  of  Ratios  of  SETA -2  Density  Data  to  J71  Model 
Plotted  as  a  Function  of  Geographic  Latitude 


3.5  Effect  of  Digital  Filtering  on  Group  II  Data 

A  quantitative  description  of  the  change  in  the  standard  deviation  of  Group  II 
data  due  to  latitude-dependent  effects  can  be  achieved  if  the  data  are  filtered  such 
that  essentially  all  wave  structures  are  passed  through  the  filter.  Preliminary 
processing  techniques  were  developed  in  an  attempt  to  reduce  ten  selected  days  of 
SETA-2' data  with  a  M15-20"-sec  filter.  Five  magnetically  quiet  days,  one  low 
activity  day  and  four  highly  disturbed  days  were  processed. 

Figure  13  compares  standard  deviations  vs  latitude  for  three  very  quiet 
(Kp  <  1)  days,  22-24  May  82,  using  both  the  "15-20"-sec  and  "50-100"-sec  filters. 
Results  are  given  for  the  MSIS77  and  J7  1  models.  The  "15-20" -sec  filter  gives 
standard  deviations  about  0.  5  ±  0.  5  percent  greater  than  the  "50 -100" -sec  filter. 

Some,  possibly  all,  of  this  very  small  difference  is  attributable  to  additional  vehicle 
dynamic  accelerations  passed  by  the  "15-20"-sec  filter.  For  latitude  bins  30*  N-60*  N 
on  the  nightside  no  data  are  plotted  for  the  "15-20"-sec  filter.  This  la  because 
noise  spikes,  removed  by  the  "50-100"-sec  filter,  were  not  rejected  by  the 
"  15 -20" -sec  filter  and  produced  large  anomalies  In  the  standard  deviations.  The 
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Figure  11c.  Standard  Deviations  of  Ratios  of  S3-4  Density  Data  to  J7  1  Model 
Plotted  as  a  Function  of  Geographic  Latitude 


present  study  indicates  that,  during  extremely  low  geomagnetic  conditions,  wave 
structures  with  periods  less  than  about  50  sec  have  very  little  effect  on  the  standard 
deviation.  Figure  13  is  consistent  with  the  results  for  the  low  Kp  group  of  Figure  11, 
in  that  it  shows  the  maximum  standard  deviation  occurs  at  high  latitudes.  However, 
for  the  lower  average  geomagnetic  conditions  (Figure  13),  the  peak  standard  devia¬ 
tion  is  about  2  percent  lower. 

Results  from  three  other  periods  studied  with  the  "15-20 "-sec  filter  are  pre¬ 
sented  in  Figure  14  in  which  the  latitudinal  distribution  of  the  standard  deviation  is 
related  to  MSIS77.  The  data  represent  12-hr  periods  each  for  low  geomagnetic 
conditions  on  10  July  82  (Kp  -  2+)  and  highly  disturbed  periods  of  12  July  82 
(Kp  ~  5+)  and  14  July  82  (Kp  ~  7+).  These  results  are  considered  tentative  because 
cf  the  difficulty  in  determining  bias  during  storm  periods  without  the  availability 

4 

of  quiet-time  data  near  these  periods,  and  the  need  to  obtain  preliminary  results 
on  a  short-term  basis.  Each  curve  is  based  on  data  from  eight  consecutive  orbits. 
For  10  July  82,  the  standard  deviation  is  less  than  5  percent  at  latitudes  below  50°  N 
on  the  dayside.  It  reaches  a  peak  of  14  percent  near  the  pole.  For  12  July  82  the 
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Figure  12,  Standard  Deviations  of  Ratios  of  S3-1  Density  Data  to  J71  Model 
Plotted  as  a  Function  of  Geographic  Latitude 


maximum  standard  deviation,  16  percent,  is  at  70-80°  N  on  the  dayside.  The  maxi¬ 
mum  for  14  July  82  is  19  percent  and  occurs  in  the  nightside  90-80°  latitude  bin. 

The  standard  deviations  exceed  that  of  10  July  at  all  latitudes  except  below  about  40° 
on  the  nightside.  Although  these  results  are  preliminary,  they  indicate  trends  simi¬ 
lar  to  those  of  Figure  11.  The  greater  high  latitude  variations  are  assumed  to  be 
related  to  the  extremely  high  level  of  geomagnetic  activity.  However,  a  much  more 
extensive  study  of  the  recently  acquired  SETA-2  data  and  the  required  processing 
techniques  using  the  M15-20"-sec  filter,  will  be  carried  out.  Such  an  effort  would 
provide,  for  the  first  time,  an  accurate  assessment  of  the  influence  of  wave 
structures  on  density  variability  over  the  range  of  conditions  encountered  by  the 
Group  II  data. 
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Figure  13.  Effect  of  Digital  Filter  Bandwidth  on  SETA-2 
Standard  Deviations  for  (a)  MSIS77  Model  and  (b)  J71  Model 
for  Quiet  Geomagnetic  Conditions 


4.  ATMOSPHERIC  VARIABILITY 

Results  given  in  Section  3  show  that  the  best  models  determine  the  average 
values  of  density  and  standard  deviation  ( lo)  to  about  ±  10  percent  and  ±  15  percent, 
respectively.  Development  of  more  accurate  models  will  require  significant  im¬ 
provements  in  our  knowledge  of  the  magnitude,  distribution  and  transport  of 
thermospheric  heating.  This  section  surveys  the  problems  involved  in  obtaining 
the  information  and  in  relating  it  to  density  variability. 

The  complex  phenomena  contributing  to  atmospheric  heating  and  dynamics  are 
summarized  in  Figure  15,  Solar  EUV,  the  major  thermospheric  heat  source,  is 
deposited  mainly  at  low  latitudes  in  the  summer  hemisphere.  The  circulation  and 
structure  of  the  thermosphere  at  low  and  middle  latitudes  are  controlled  primarily 
bv  this  heating.  Upward  propagating  gravity  waves  and  tides  provide  another  energy 
source  effective  mainly  at  low  and  middle  latitudes.  At  auroral  latitudes  electro¬ 
dynamic  coupling  occurs  between  the  magnetosphere,  ionosphere  and  neutral 
atmosphere.  Both  the  magnitude  and  spatial  extent  of  the  high- latitude  heating  are 
extremely  variable  and  are  related  to  the  level  of  geomagnetic  activity.  Energetic 
charged-particle  precipitation  from  the  magnetosphere  and  Joule  (resistance'  heating 
due  to  electrical  current  systems  within  the  ionosphere  are  important  sources 
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of  energy.  These  localized  heat  sources  can  cause  large  and  sudden  density 
variations,  perturb  circulation  patterns  and  launch  gravity  waves.  Also, 
ionospheric  plasma,  driven  by  the  magnetospheric  convection  electric  field, 
transfers  momentum  to  the  neutral  particles  and  can  produce  strong  winds.  The 
wind  systems  induced  by  high  latitude  processes  can  significantly  alter  thermospheric 
density,  composition  and  temperature  on  a  global  scale  during  disturbed  geomag¬ 
netic  conditions. 
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Figure  14.  Standard  Deviations  of  Ratios  of  SETA-2  Data  to  MSIS77  IVlodel  for 
Low,  Moderate  and  High  Geomagnetic  Periods  Plotted  as  a  Function  of  Geo¬ 
graphic  Latitude 
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Figure  15.  Schematic  Representation  of  Heat  Sources  for  the  Upper  Atmosphere 


Density  in  the  lower  thermosphere  is  known  to  vary  with  the  following 
parameters: 

(1)  Solar  flux  (solar  cycle  and  dally  components) 

(2)  Geomagnetic  activity 

(3)  Local  time 

(4)  Day  of  year 

(5)  Latitude 

(G)  Longitude 

(7)  Wave  structures. 

Problems  in  achieving  improved  modelling  of  each  of  these  factors  are  dis¬ 
cussed  below.  However,  it  must  be  understood  that  the  energetics,  dynamics, 
chemistry,  composition  and  resulting  structure  of  the  atmosphere  are  all  functions 
of  these  factors  and  are  so  intricately  related  that  it  is  not  really  possible  to  discuss 
each  in  isolation. 

Solar  Flux  -  The  ^  index  used  in  all  current  models  does  not  necessarily 
represent  the  complex  mechanisms  of  interactions  between  the  EUV  flux  and  the 
thermosphere.  Since  it  has  been  found  to  generally  reflect  variations  in  the 
thermospheric  energy  input,  it  is  routinely  used  as  a  readily  available,  but  imper¬ 
fect,  indicator  of  solar  activity.  It  presents  an  additional  problem  for  near  real 
time  operations.  The  solar-cycle  dependent  value  of  FlQ  ^  used  in  models  is  based 
on  data  some  three  months  either  side  of  the  time  of  interest.  Studies  to  determine 

12 

the  relationship  between  F^  ^  and  satellite  measurements  of  EUV  are  in  progress. 

An  EUV  reference  spectrum  for  1976-1981  has  been  developed.  However,  extrapola¬ 
tion  outside  this  period  is  not  recommended.  Correlations  have  been  found  to  be 
non-linear  and  both  wavelength  and  solar-cycle  dependent.  In  order  to  achieve  the 
goal  of  ultimately  utilizing  direct  measurements  of  the  solar  EUV  flux  in  thermos¬ 
pheric  models,  an  extensive  measurement  and  data  analysis  program  will  be  re¬ 
quired.  Methods  of  modelling  UV  irradiance  based  on  detailed  ground-based  ob¬ 
servations  of  the  distribution  of  specific  active  regions  across  the  solar-disk  are 

13 

also  being  investigated.  Some  promising  results  have  been  obtained. but  these 
studies  are  in  preliminary  stages. 

Geomagnetic  Activity  -  Large  (up  to  ~  100  percent)  and  short-term  perturbations 
in  neutral  density  and  winds  (~  1  km/ sec)®'  14  have  been  observed  at  high  latitudes 
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in  response  to  geomagnetic  storms.  Both  the  amplitude  and  phase  of  these  fluctua¬ 
tions  must  be  accurately  modelled  to  significantly  Improve  operational  predictive 
capabilities.  The  problem  is  complicated  because  the  storms  are  caused  by  high 
latitude  heating  processes  that  are  not  well  understood  and  the  density  response 
depends  on  latitude,  longitude,  season,  local  time  and  solar  cycle  as  well  as  storm 
intensity, 

The. planetary  3 -hr  Kp  or  ap  index  used  as  an  indicator  in  most  models  does 
not  necessarily  represent  the  physical  mechanisms  responsible  for  density  varia¬ 
tions.  However,  It  is  the  only  routinely  available  Index  for  geomagnetic  activity. 
Based  on  a  network  of  midlatitude  stations,  it  is  particularly  limited  when  des¬ 
cribing  data  at  high  latitudes. 

Correlations  of  geomagnetic  activity  with  the  interplanetary  magnetic  field  (B^) 

magnitude  and  direction  and  the  solar  wind  velocity  (V)  are  being  studied.  Relations 

of  the  form  B.  V  where  0.  85  5  m  =  2  and  1  >  n  5  2  have  been  shown  to  provide  a 

measure  of  energy  transferred  from  the  solar  wind  to  the  magnetosphere.  These 

parameters  also  correlate  highly  with  the  auroral  zone  electrojet  indices  AE  and 
15 

AL.  Eventually  it  may  be  possible  to  utilize  an  index  derived  from  these  solar 

wind  and  interplanetary  magnetic  field  measurements  as  a  precursor  of  geomagnetic 

activity.  Considerable  research  to  develop  improved  formulations  relating  any  new 

index  to  the  complicated  density  response  would  still  be  required. 

1 6 

The  Dynamics  Explorer  satellites  are  monitoring  relevant  parameters  in¬ 
volved  in  coupling  between  the  magnetosphere,  ionosphere  and  thermosphere.  These 

measurements  are  being  obtained  at  altitudes  above  300  km.  They  provide  direct 

17 

inputs  for  evaluation  and  improvement  of  dynamic  global  circulation  models  and 
should  lead  to  considerable  progress  in  understanding  the  physics  of  geomagnetic 
storms.  However,  the  resultant  improvement  to  be  realized  in  the  specification  of 
both  the  amplitude  and  phase  of  density  variations  due  to  a  particular  storm  cannot 
be  estimated  at  present.  Also,  the  extrapolation  of  Dynamics  Explorer  data  to  alti¬ 
tudes  below  200  km  provides  significant  difficulties.  There  are  few  measurements 
of  energetics  and  dynamics  in  this  region  to  guide  theoretical  developments. 

Local  Time  -  The  dominance  of  the  semidiurnal  tide  below  200  km  has  recently 

18  19  20 

been  revealed  by  low-altitude  satellite  density  measurements.  At  present, 

the  referenced  Jacchla  models  do  not  include  this  contribution  to  local  time  varia¬ 
tions.  Extensive  theoretical  progress  has  been  made  in  describing  both  the  diurnal 

and  semidiurnal  components  as  a  function  of  altitude,  latitude,  season  and  solar 
21  22 

flux.  '  In  general,  these  results  show  good  agreement  with  experimental  data 
above  150  km.  At  lower  altitudes  significant  day-to-day  irregularities  persist. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  47.) 


37 


This  Is  related  to  the  Interaction  betwoen  local  winds  and  tidal  processes.  These 
complicated  dynamics  make  It  difficult  to  establish  the  average  tidal  behavior. 

Over  moat  of  the  region  below  140  km  the  largest  density  perturbations  generally 
occur  at  equatorial  latitudes  for  the  diurnal  and  semidiurnal  tides.  At  140  km  the 
density  uncertainty,  depending  on  the  phase  of  these  two  components,  can  be  between 
15  percent  and  40  percent  (J.  M.  Forbes,  private  communication). 

Day  of  Yoar  -  Density  variations  with  annual  and  sub-annual  periodicities,  as 
distinguished  from  periodic  changes  connected  with  the  earth  and  sun.  are  a  regular 
feature  of  the  upper  atmosphere.  The  "semi-annual  variation"  Is  generally  charac¬ 
terized  In  current  models  by  maxima  near  equinox  and  minima  near  solstice.  At 

the  present  time,  the  causo  of  the  semi-annual  variation  has  not  been  unambiguously 

23 

resolvod.  Ching  und  Chiu  have  shown  that  EUV  heating  has  the  proper  character¬ 
istics  to  be  important  at  altitudes  above  150  km.  Mayr  and  Volland2*  have  shown 
that  sources  in  addition  to  EUV  heating  are  required  to  account  fully  for  this  varia¬ 
tion.  High  latitude  heat  sources,  including  auroral  heating  related  to  the  semi¬ 
annual  component  of  the  occurrence  of  magnetic  storms  have  been  Invoked.  Walter- 
schcid25  recently  has  considered  the  semi-annual  density  variation  as  a  conduction 
mode  of  oscillation  essentially  driven  by  Joule  heating  at  high  latitudes.  This 
theoretical  analysis  shows  good  agreement  with  amplitude  and  phase  vs  altitude 
and  latitude  predicted  by  MSIS.  Additional  studies  are  needed  to  attain  a  more 
complete  understanding  of  the  semi-annual  effect.  MSIS  (and  other  models)  show 
a  larger  amplitude  maximum  in  October  than  in  April.  Measurements  using  the 
satellite  orbital  decay  technique2®  indicate  that  during  low  solar  flux  conditions  the 
reverse  is  true. 

Latitude  -  In  addition  to'the  latitude  dependences  associated  with  each  of  the 
above  phenomena,  present  models  also  recognize  seasonal-latitudinal  variations. 
Jacchln7  has  described  contributions  to  this  effect  at  thermospheric  heights  due  to 
both  the  seasonal-latitudinal  component  of  the  diurnal  temperature  variation  (strongly 
height  dependent)  and  an  intrinsic  seasonal-latitudinal  variation  (not  height  dependent) 
with  origin  at  the  lower  boundary  of  the  thermosphere.  In  addition,  below  about 
150  km,  there  Is  a  variation  due  to  the  stratospheric  and  mesospheric  seasonal 
variation  which  "spills  over"  into  the  lower  thermosphere. 


23.  Ching,  13.  K. ,  and  Chiu,  Y,  T.  (1072)  Annual  and  sub-annual  effects  of  EUV 

heating-I.  Harmonic  analysis,  Planet.  Space  Scl.  20j  1745. 

24.  Muyr,  II.  G. ,  und  Volland,  II.  (1972)  Theoretical  model  for  the  latitude  depend¬ 

ence  of  the  annual  and  semiannual  variations,  J.  Oeophys.  Res.  ,JXi0774, 

25.  Walterscheid,  R.  L.  (1982)  Semiannual  Oscillation  in  the  Thermosphere  as  a 

Conduction  Mode,  Aerospace  Corp.  ftept.  No.  ATR-b2  (0555)  -  I. 

20,  Walker,  D,  M.  C.  ( 1078)  Variations  In  air  density  from  January  1972  to  April 
1075  at  heights  near  200  km,  Planet.  Space  Scl.  Jffit 20 1 . 
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Longitude  -  Satellite  composition  measurements  have  provided  evidence  of  the 
Influence  of  dynamics  on  longitudinal  variations.  Density  maxima  for  heavy  con¬ 
stituents  occur  at  the  longitude  of  the  north  and  south  magnetic  poles.  The  behavior 
of  the  lighter  constituents  is  negativoly  correlated.  This  behavior  is  related  to 
plasma  density  and  motion  influenced  by  earth's  magnetic  field.  Ion-neutral  coupling 
is  therefore  Involved  In  determining  the  longitudinal  distribution  of  the  neutral  den¬ 
sity.  The  MSIS  1 07 9®  model  included  terms  dependent  upon  geographic  latitude, 
longitude,  universal  time,  solar  flux  and  geomagnetic  activity  to  account  for  this 
variation.  However,  the  standard  deviation  of  the  data  shown  in  Table  &  waa  about 
0.  3  percent  higher  with  this  model  than  with  the  MSIS  1079  model  which  waa 
longitudinally  averaged,  indicating  the  difficulty  in  modelling  this  variation. 

Waves  -  Large-scale  gravity  waves  having  periods  different  from  those  of  tidal 
waves  and  not  directly  related  to  the  earth- sun- moon  geometry,  also  occur  in  the 
atmosphere.  These  may  be  generated  either  locally  in  the  thermosphere  or  at 

lower  altitudes.  Gravity  waves  with  the  largest  amplitudes  observed  at  high  latl- 

27 

tudcs  were  detected  by  the  accelerometers  on  the  OV1-15  satellite.  These  were 
related  to  the  intensity  of  high-latitude  heating  as  indicated  by  the  auroral  electro- 

ofl 

jet  index.  The  Kp  index  was  found  to  be  a  poor  indicator  of  wave  generation. 

Evidence  of  wave  structure  with  location  dependent  on  earth's  tropospheric  weather 

29 

patterns  has  been  found  in  the  analysis  of  Atmosphere  Explorer  ion  gauge  data. 

A  comprehensive  survey  of  wave  structures  in  neutral  constituents  was  made  using 

30 

330  orbits  of  data  from  the  open-source  neutral-mass  spectrometer  on  AE-C. 

The  highest  wave  amplitudes  (up  to  55  percent  peak-to-peak)  and  greatest  number 
of  occurrences  were  in  both  the  northern  and  southern  high-latitude  regions.  These 
were  associated  with  a  high-latitude  conjugate  phenomenon  such  os  Joule  heating 
and/or  particle  precipitation.  At  ±  25*  magnetic  latitude,  secondary  maxima  in 
low  amplitude  waves  (5-15  percent)  were  attributed  to  a  near  equatorial  source. 

These  wave  phenomena  are  not  incorporated  into  present  models.  While  some 
of  the  features  of  wave  structures  associated  with  geomagnetic  activity  are  becoming 
known,  a  capability  for  accurately  predicting  such  behavior  may  not  be  possible  in 
the  foreseeable  future.  Also,  planetary  and  gravity  waves  with  sources  in  the 
troposphere  or  stratosphere  appear  to  propagate  into  the  thermosphere  and  dissipate 


27.  Marcos,  F.  A, ,  and  Champion,  K.S,  W.  (1972)  Gravity  waves  observed  in  high 

latitude  neutral  density  profiles,  Spacs  Res.  £&79l. 

28,  Forbes,  J.M.,  and  Marcos,  F.A,  (1973)  Thermospheric  density  variations 

associated  with  auroral  eloctrojet  activity,  J,  Geophvs.  Res.  |£|384l. 

20,  Rice,  C.  J. ,  and  Sharp,  L,  R.  (1077)  Neutral  atmospheric  waves  in  the 

thermosphere  and  tropospheric  weather  systems,  Oeophys.  Res,  Lett.  ^3 15a. 

30,  Potter,  W.  E.,  Kayser,  D.  C, ,  and  Mauereberger,  K.  (1978)  Direct  measure- 
.  ments  of  neutral  wave  characteristics  in  the  thermosphere,  J.  Geophvs,  Res, 
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energy  by  viscous  interactions.  Their  influence  on  heating  and  transport  must  be 
investigated. 

Table  8  summarizes  the  approximate  magnitude  and  time  scale  for  density 
increases  associated  with  each  of  the  above  phenomena.  Results  are  given  for  an 
altitude  of  200  1cm.  Magnitudes  are  representative  of  those  in  current  models 
except  for  the  wave  value,  which  is  based  on  measured  data. 


Table  8.  Thermospheric  Density  Variability  at  200  km 


Effect 

Ap 

200  km 

Time 

Scale 

1. 

Flux  (Solar  Cycle) 

100% 

Years 

2. 

Flux  (Daily) 

10% 

Day 

3. 

Geomagnetic  Activity 

100% 

Hours 

4'. 

Local  Time 

15% 

Hours 

5. 

Semi-Annual 

40% 

Months 

6. 

Latitude 

25% 

Months 

7. 

Longitude 

5% 

Day 

8. 

Waves 

50% 

Minutes 

5.  RECOMMENDATIONS 

'C 

Sections  3  and  4  have  shown  that  the  existing  need  for  very  accurate  knowledge 
of  both  absolute  density  and  its  variability  require  a  long-term  program  involving: 

(a)  More  realistic  indices  to  describe  energy  inputs  related  to  solar  EUV 
and  auroral  processes. 

(b)  An  extensive  measurements  program  to  provide  a  basis  for  developing 
accurate  relationships  between  energy  inputs  and  density  variability.  N 

{ c)  Development  of  accurate  dynamic  global  circulation  models  that  can  better 
determine  how  the  thermosphere  is  affected  by  internal  and  external  influences. 

Many  of  the  efforts  directed  toward  replacing  F1()  7  and  Kp  were  discussed  at 
the  Satellite  Drag  Workshop.  31  It  was  concluded  that  more  realistic  solar  and  geo¬ 
magnetic  indicators  are  still  at  least  several  years  away.  Until  better  indices 
become  available,  there  will  be  a  fundamental  limitation  on  model  accuracies. 


31.  Space  Environment  Services  Center  (1982)  Proceedings  of  a  Workshop  on 
Satellite  Drag,  March  18-19,  1982,  U.  S.  Dept,  of  Commerce, 

Boulder,  Colorado. 


40 


There  remains  u  critical  need  for  data  at  all  altitudes  In  the  lower  thermosphere 
The  problem  is  particularly  acute  at  high  latitudes  and  below  150  km.  Atmosphere 
Explorer  data  were  obtained  mainly  above  150  km.  Also.  AE-D,  the  polar  orbit 

mission,  lasted  only  a  few  months  due  to  a  spacecraft  power  failure.  The  Space 
32 

Science  Board  has  pointed  out  the  importance  of  determining  the  mean  circulation 

und  temperature  structure  as  well  as  changes  in  the  tidal  structure  and  changes  in 

geomagnetic  activity  to  understand  the  region  below  150  km. 

Global  circulation  models  are  being  continually  updated  and  Improved.  A 

recent  revision  Includes  the  incorporation  of  the  effects  of  the  high-latitude  convec- 
33 

tion  electric  field.  Dynamics  Explorer  data  will  permit  significant  advances  in 
magnetic  disturbance  morphology  at  altitudes  above  300  km.  However,  it  should  be 
at  least  another  5-10  years  before  these  dynamic  models  can  achieve  the  accuracies 
of  present  empirical  models.  Additional  extensive  effort  beyond  this  time  frame 
will  be  needed  to  extend  these  models  accurately  at  altitudes  below  200  km  and  to 
demonstrate  better  accuracy  than  empirical  models. 

It  appears  that,  at  least  through  this  decade,  there  will  not  be  an  alternative 
to  direct  measurements  for  knowledge  of  absolute  density  and  Its  variations.  Re¬ 
mote  sensing  techniques  for  global  monitoring  are  being  investigated,  but  even  if 
implemented  successfully,  they  may  not  achieve  the  required  accuracies.  Even¬ 
tually,  as  models  and,  possibly,  remote  sensing  techniques  improve,  directly 
measured  data  may  be  needed  only  for  normalization,  updating  and  fine-tuning  of 
model  predictions. 

The  remainder  of  this  section  deals  with  three  requirements  to  improve  under¬ 
standing  of  absolute  density  and  its  variability:  first,  a  technique  to  provide 
required-accuracy  density  measurements  using  the  satellite  accelerometer  experi¬ 
ment;  second,  utilizing  this  technique  an  extensive  measurements  program  to 
provide  the  basis  for  development  of  more  accurate  models;  and  third,  other 
needed  supporting  efforts. 

S.l  Absolute  Density  Determination 

An  improved  method  of  determining  absolute  density  can  be  obtained  by 
exploiting  accelerometer  capability  to  provide  extremely  accurate  drag  measure¬ 
ments.  Density  is  extracted  from  the  drag  data  by  the  relation: 


32.  Committee  on  Solar- Terrestrial  Research  (1981)  Solar- Terrestrial  Research 

for  the  1980's,  National  Academy  Press,  Washington,  D. C. 

33.  Roble,  R.  G. ,  Dickinson,  R.  E. ,  and  Ridley,  E.  C.  (1982)  Global  circulation 

and  temperature  structure  of  thermosphere  with  high  latitude  plasma 
convection,  J.  Geophys,  Res.  J&tl599. 
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(3) 


2Muq 

p  s  - - — - 

A  C^v  cob  0 

where  p  Is  density,  la  the  measured  drag  acceleration,  0  is  the  angle  between  ? 
the  accelerometer  and  drag  vector  and  M,  A,  and  V  are,  respectively,  the 
satellite's  mass,  cross  sectional  area,  drag  coefficient  and  the  atmospheric  mass 
velocity  (relative  to  the  satellite).  The  required  experiment,  which  could  routinely 
measure  absolute  density  to  about  ±  5  percent,  should  have  the  following  charac¬ 
teristics: 

(1)  A  spherical  shape  for  accurate  area  determination  and  improved 
theoretical  drag  coefficient  estimate. 

(2)  An  elliptical  orbit  to  permit  calculation  of  accelerometer  noise  and 
bias  near  apogee. 

(3)  An  accurate  attitude  system  to  reduce  errors  in  0. 

(4)  A  minimum  of  extraneous  accelerations  due  to  vehicle  dynamics. 

(5)  A  lifetime  of  several  months  in  a  relatively  low  orbit  (below  150  km) 
to  acquire  data  over  a  wide  range  of  variables. 

(6)  Incorporation  of  a  wind  sensor  to  reduce  errors  in  determining  the 
atmospheric  velocity  relative  to  the  satellite. 

(7)  Incorporation  of  a  calorimeter  experiment  to  further  improve  drag 
coefficient  estimates  (A.  Prag,  private  communication). 

A  mass  spectrometer  should  also  be  a  part  of  the  above  experiment.  Although 
it  may  not  provide  as  accurate  total  density  data  as  the  accelerometer  experiment, 
the  composition  data  are  invaluable  for  interpreting  variations  in  density.  These 
experiments  also  complement  each  other  since  the  accelerometer  can  provide 
vertical  and  cross-track  winds  and  the  mass  spectrometer  can  provide  along-track 
winds. 

5.2  Measurement!  Program! 

The  satellite  accelerometer  experiment  is  crucial  for  advancing  our  knowledge 
of  density  and  its  variability  in  the  lower  thermosphere.  To  date  it  has  provided 
the  only  comprehensive  set  of  measurements  during  high  solar  flux  conditions  as 
well  as  an  extensive  low  solar  flux  data  base.  As  shown  in  Sections  3  and  4,  addi¬ 
tional  extensive  measurements  are  needed  to  permit  continued  improvements  in 

32 

operational  capabilities.  Also,  the  Space  Science  Board  has  stressed  the  need 
for  continued  spacecraft  Investigations  of  the  upper  atmosphere  throughout  the 
1980's  to  ensure  adequate  geographic,  local  time,  seasonal  and  solar-cycle  coverage. 

The  first  set  of  measurements  suggested  is  an  extension  of  the  Group  II  data 
set  to  lower  solar  flux  conditions  as  follows: 
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(1)  An  accelerometer  flown  during  1983  has  provided  data  during  intermediate 
solar  flux  conditions.  These  data  must  be  incorporated  into  the  density  data  base 
and  analyzed. 

(2)  An  opportunity  exists  to  fly  an  accelerometer  on  a  host  vehicle  during  1985. 

A  partially  fabricated  instrument  could  be  completed  in  time  to  meet  the  flight 
schedule.  Correlative  measurements,  particularly  a  mass  spectrometer,  should 
also  be  flown  on  this  mission.  Other  desired  payloads  include  a  wind  sensor  and 

a  calorimeter.  This  mission  would  provide  data  near  solar  minimum  and  would 
furnish  systematic  measurements  of  high-latitude  dynamic  processes  that  were  not 
obtained  with  AE-D. 

These  two  flights,  combined  with  the  Group  II  data,  would  provide  a  consistent 
data  set  to  improve  formulations  of  density  variability  as  a  function  of  solar  flux, 
geomagnetic  activity,  day  of  year,  longitude  and  wave  structures  as  a  function  of 
latitude  above  170  km.  This  program  should  be  extended,  using  the  payload  sug¬ 
gested  for  item  2  above,  with  heat  source  measurements,  to  obtain  year-round  data 
in  the  northern  hemisphere  in  (1)  a  noon-midnight  orbit  and  (2)  a  dawn-dusk  orbit. 
These  missions  should  overlap  by  a  few  months  to  permit  more  detailed  analysis  of 
local  time  variations  and  geomagnetic  storm -related  density  and  circulation  varia¬ 
tions. 

The  next  set  of  measurements  should  provide  an  expanded  study  of  the  dynamics 
of  the  lower  thermosphere.  The  absolute  density  experiment  of  Section  5.  1  with 
composition  and  wind  data  should  routinely  explore  density  structure  ataltitudes 
as  Low  as  130  km.  A  pair  of  simultaneous  satellites  with  lifetimes  greater  than 
one,,year  is  suggested.  The  spacecraft  would  consist  of  a  low-  and  a  high-inclination 
mission.  The  concept  is  similar  to  the  overlapping  AE-D  and  AE-E  missions  con¬ 
cept  which  was  carried  out  for  only  a  short  time.  The  objective  would  be  to  pro¬ 
vide  accurate  density  data  at  critical  lower  altitudes.  Such  missions  will  require 
trade-offs  between  available  propulsion  capability,  perigee  and  apogee  heights  and 
experiment  operations.  Follow-on  missions  are  needed  for  solar-cycle  coverage. 
These  missions  should  also  be  coordinated  with  other  onboard  experiments  or 
satellite  and  ground-based  experiments  that  monitor  solar  energy  inputs,  ionospheric 
plasma,  electric  fields  and  charged  particle  precipitation.  Tethered  satellite 
missions  should  also  be  considered,  when  available. 

5.3  Supporting  Studies 

Other  efforts  that  can  lead  to  an  improved  capability  to  specify  the  density  at 
satellite  altitudes  include  the  following: 


43 


(1)  Below  about  250  km  there  ia  very  little  data  on  atmoapheric  dynamics. 
Detailed  studies  of  the  morphology  of  density  and  wind  during  storm  periods 
encountered  by  SETA-1  are  being  carried  out.  These  studies  will  also  incorporate 
temperatures  derived  from  incoherent  scatter  radar  data.  The  above  efforts  form 
the  basis  for  development  of  the  required  more  extensive  analytes.  These  analysis 
efforts  involving  valuable  data  in  regions  of  particular  interest  should  be  continued 
with  the  more  recent  data  sets. 

(2)  New  indices  for  solar  and  geomagnetic  activity  should  be  tested,  when 
available,  using  the  present  Group  I  and  Group  II  accelerometer  data  sets  (assuming 
the  indices  are  based  on  data  available  during  the  time  periods  of  the  density 
measurements). 

(3)  Both  empirical  and  dynamic  models  should  be  improved  using  results  from 

above  analyses.  . 

(4)  The  feasibility  of  using  auroral  imaging  techniques  or  heating  measure¬ 
ments  to  relate  auroral  zone  heating  to  density  variability  should  be  investigated. 

In  particular,  it  should  be  determined  whether  regions  of  intense  heating  can  be 
related  to  wave  generation. 


6.  CONCLUSIONS 

......  .A  ) 

Satellite  density  measurements  using  accelerometers  have  provided  an  extensive 
description  of  the  variability  of  the  neutral  atmosphere.  Mean  densities  are 
believed  to  be  given  within  ±  10  percent  by  the  most  accurate  models.  This 
uncertainty  can  be  reduced  to  about  ±  5  percent  by  an  accelerometer  experiment 
such  as  described  in  Section  5.  1.  Low  solar  flux  results  derived  from  AE-C,  -D  * 
and  -E  and  S3-1  satellites  indicate  that  the  best  models  have  a  standard  deviation 
of  about  ±  15  percent.  For  the  higher  solar  flux  data  from  SETA-1,  SETA -2  and 
S3 -4  the  standard  deviation  was  near  9  percent.  This  lower  value  is  attributed  to 
the  additional  smoothing  due  to  the  50-100  sec  filter  and  the  high  correlation  of 
altitude,  latitude  and  local  time  for  this  data  set.  The  measurements  did  show  an 
increased  model  uncertainty  at  high  latitudes. 

The  following  efforts  can  be  easily  accommodated  into  existing  programs  and 
-provide  important  results  needed  to  augment  and  define  future  activities: 

(1)  A  consistent  set  of  data  with  excellent  solar-cycle  coverage  can  be  achieved 
by  a  1985  flight  of  an  accelerometer  on  a  host  vehicle.  The  results  would  extend 
the  present  Group  II  data  base  to  low  solar  flux  conditions  and  provide  a  unique 
basis  for  evaluation  and  improvement  of  density  models. 
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(2)  Density  and  wind  measurements  from  the  SETA-1  and  -2  experiments 
already  provide  a  unique  opportunity  to  obtain  fundamental  new  results.  Data 
analysis  efforts  are  being  made  to  understand  the  dynamic  processes  influencing 
density  variability.  Such  studies  will  lead  to  improved  modelling  of  density  in  the 
thermosphere.  Additional  data  analysis  support  would  produce  these  results  on 
a  timely  basis. 

The  accuracy  that  models  can  ultimately  achieve  cannot  be  specified  at  present. 
Further,  sufficient  data  are  not  yet  available  to  fully  assess  present  model  errors. 
The  problems  involved  in  improving  models  and  an  indication  of  some  of  the  on¬ 
going  efforts  directed  toward  their  resolution  were  summarized  in  Section  4. 

Major  problem  areas  are  wave  structures  and  high-latitude  processes  involved 
during  geomagnetic  storms.  Significant  reductions  in  density  uncertainties  appear 
to  be  achievable  but  require  long-term  programs  involving  synoptic  measurements 
coordinated  with  theoretical  studies  and  modelling  efforts. 
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